Dislocation induced anomalous softening of solid helium 
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The classical motion of gliding dislocation lines in slip planes of crystalline solid helium leads to 
plastic deformation even at temperatures far below the Debye temperature and can affect elastic 
properties. In this work we argue that the gliding of dislocations and plasticity may be the origin 
of many observed elastic anomalies in solid "^He, which have been argued to be connected to super- 
solidity. We present a dislocation motion model that describes the stress-strain r-e curves and work 
hardening rate dr/de of a shear experiment performed at constant strain rate e in solid helium. The 
calculated dr/de exhibits strong softening with increasing temperature due to the motion of dislo- 
cations, which mimics anomalous softening of the elastic shear modulus fi. In the same temperature 
region the motion of dislocations causes dissipation with a prominent peak. 

PACS numbers: 61.72.Hh, 67.80.B-, 67.80.bd 



The observation of softening of the low-temperature 
shear modulus in solid ^He with increasing temperature 
around 100 mK has been taken as evidence for anoma- 
lous elastic properties tied to supersolidity [J-Q. Mea- 
surements of the shear modulus and the resonant period 
of torsional oscillator of solid helium showed no period 
drop in the case of "^He, while both ^He and '*He systems 
showed softening of the shear modulus in the same tem- 
perature range. This result suggested the importance of 
the role of quantum statistics and the role of defects and 
moving dislocations on supersolidity [HI . The connection 
of the motion of dislocations with supersolidity is not cer- 
tain, whereas the connection of the glide of dislocations 
with plastic deformation of the crystal is well established 
0. The theory of strain hardening in materials science 
predicts the stress-strain r-e curve and allows the study 
of dislocation motion and plasticity 0, II] • Such disloca- 
tion motion would strongly affect the determination of 
elastic properties in solid helium, a crystal with a very 
soft elastic shear modulus 

In this Letter, we present a new perspective on the 
origin of reported anomalous elastic properties of solid 
helium. The glide of dislocations in a dislocation forest 
is well understood for typical metals and can lead to large 
changes of order 60% in the work hardening rate (WHR) 
[lO|. Here, we study the dynamics of classical dislocation 
motion and plastic behavior of solid helium. At lower 
temperatures one might expect that quantum dislocation 
motion will play a significant role [lllll3- In the case of a 
DC shear strain rate experiment, e = const., with shear 
stress T, we predict the softening of the WHR, dr/de, 
with increasing temperature T due to dislocation glide 
and the creation and multiplication of dislocation loops. 
In the limit of e — > the WHR approaches the elastic 
shear modulus dr/de — > /i and exhibits purely elastic be- 
havior. Our dislocation model also predicts that higher 
■^He impurity concentration in solid ^He pushes the on- 
set of the magnitude change of dr/de to higher temper- 



atures, because more '^He atoms pin more dislocations, 
and the pinning increases the number of immobilized dis- 
locations, causing the crystal to harden. At the same 
time, the zero-temperature value of the WHR is unaf- 
fected by ^He atoms, whereas its high-temperature value 
is reduced by a lower amount with higher '^He concentra- 
tions. Finally, the WHR decreases with increasing strain 
at finite temperature and attains the elastic shear mod- 
ulus value at absolute zero temperature. 

Our proposed scenario for dislocation dynamics in solid 
^He differs from the prevailing view of boiling off of '^He 
impurities from dislocation lines with increasing T, which 
act as pinning sites, and the subsequent vibration of 
longer dislocation segments pinned only by the disloca- 
tion network main result of this work is the 
prediction of a dissipation peak caused by plastic defor- 
mation of the solid in the same temperature region where 
dr/de changes most rapidly. All model predictions are 
remarkably similar to existing shear measurements with 
applied AC strain rate [ll-Q and will provide a stringent 
test for DC strain rate experiments of solid helium. 

We model the elastic and plastic properties of a poly- 
crystalline sample of ^He in the presence of a uniform ap- 
plied shear stress Tapp (see Fig.[IJa)). Plastic deformation 
arises when the dislocations in the solid start to glide. In 
the present model, we adopt a statistical representation 
of a group of dislocations, which neglects individual dis- 
location characteristics such as polarity, line orientation, 
etc. The dislocation network is then represented by a 
continuous distribution of dislocations and characterized 
by a linear density p (total length of dislocation lines over 
the sample volume). The total dislocation density in a 
material consists of both mobile and stored (immobile) 
dislocations 

Ptot — Pmob + P stored- (1) 

Immobile dislocations may be glissile, that is, pinned by 
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defects, or sessile, that is, lying out of their habit plane. 
For example, glissile dislocations can be pinned when dec- 
orated with ■^He atoms as sketched in Fig. [lib). Both 
Pmob and Pstored cvolve as deformation proceeds. 

During the time Ai glide increases the total dislocation 
density according to l8| 



Aptot = /SpmobvAt/L, 



(2) 



where L ~ ^/y^Ptot is the mean-free path of the mobile 
dislocations, v is their average speed, and the geometric 
coefficient < /3 < 1 is indicative of how fast the density 
increases as dislocation loops propagate The incre- 
ments in stored and mobile dislocation densities are gen- 
erally proportional to Aptot , as portions of an expanding 
dislocation loop can react with other loops and become 
permanently immobilized. Combining Eqns. ([ij and ([2|), 
we arrive at Apmob = cAptot and Ap stored = {l-c)Aptot, 
with < c < 1. 

For simplicity, we assume an average speed v for all 
mobile dislocations in the sample volume and that their 
speeds are sufficiently low, such that their motion is con- 
trolled by thermal activation obeying an Arrhenius law. 
In this regime, mobile dislocations glide from one pin- 
ning point (^He or intersecting dislocations) to another. 
Their motion is resisted by Tree, however, a combina- 
tion of thermal energy and mechanical work in excess of 
Tres can overcomc the Gibb's free energy of activation, 

AG = Ui — TeffV, with effective stress Teff = Tapp — Tres 

and activation volume V . Accordingly, the average dis- 
location speed as a function of stress and temperature is 
given by J,, M 



^biyocxp [-AG/kBT]. 



(3) 



The magnitude of the Burgers vector is 6 = |b|, the Debye 
frequency is v^, the pinning energy is Ui, and the active 
volume is V — L. We model the resistance stress by 

m 



TPe 



(4) 



with Peierls stress Tp^ieris, elastic shear modulus p, and 
typical coefficient 0.1 < a < 2 [l^. In mean- field theory, 
the second term in Q arises due to long-range interac- 
tions with other dislocations in the background. 

The total applied strain is accommodated by both elas- 
tic and plastic strain, etot — Ce + fp. The plastic strain 
rate due to glide obeys jBI 



hv p„ 



.oh- 



The corresponding applied shear stress Tapp 
Hooke's law and obeys 



Since the plastic work Wp equals Tapp tp, one should note 
that dislocation glide yields a finite plastic strain rate kp 
that gives rise to finite dissipation [Wp). 





FIG. 1: (Color online), (a) Illustration of the applied shear 
stress on the solid ^He sample with curved lines indicating 
grain boundaries; (b) magnification of circled region in (a) 
and plot of a mobile dislocation moving on the slip plane 
under external load. Black dots indicate '^He pinning points, 
solid lines are mobile dislocations, and dashed lines are stored 
(immobile) dislocations. 




FIG. 2: (Color online), (a) Stress-strain curves for difi'erent 
temperatures at applied strain rate e = 0.1 with initial 
total and mobile dislocation densities 10^ m~^ and 10* m~^, 
respectively. The filled and open circles bracket the transition 
zone between elastic and plastic behavior, (b) Plot of the cor- 
responding dislocation density-strain curves of (a). As tem- 
perature decreases, the transition point is delayed to higher 
shear strain and the transition zone narrows and sharpens. 



The stress-strain curves in Fig. (SJa) show plastic be- 
havior (deviation from the linear stress-strain curve) 
above a critical applied strain. The transition zone from 
elastic to plastic behavior is indicated by filled and open 
circles along the stress-strain curves. The applied stress 
^^•j effectively lowers the barrier given by the Gibb's free 
energy AG = Ui + TresV - TappV . When Tapp < Tc 
follows only a negligible amount of the dislocations can ther- 
mally overcome the "effective" barrier. The dislocation 
lines are essentially pinned and the sample remains elas- 
(6) tic. When Tapp > Tc an increasing amount of disloca- 



tions starts to glide, which produces plastic strain Cp for 
the same applied stress (deviation from elastic behavior). 
Accordingly dislocation loops grow in size, thus increas- 
ing the dislocation density shown in Fig.l^Jb). Tempera- 
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FIG. 3: (Color online), (a) Work hardening rate versus tem- 
perature for diflerent strains. The parameters are the same 
as in Fig. [2] (b) corresponding plastic dissipation of (a) . The 
dashed line at T = 20 mK marks the cutoff of a typical low 
temperature experiment. For thermally activated dislocation 
glide the work hardening rate approaches the elastic shear 
modulus at T = K and no plastic work is dissipated. 



ture comes into play in the elastic-plastic cross-over, since 
our model describes dislocation gliding as a thermally ac- 
tivated process in Eq. ([3]): At higher temperature, more 
dislocations can overcome the barriers. This corresponds 
to a lower critical strain and broader transition shown in 
Fig. [2ja) at high temperature. The dislocation density 
then increases at a lower applied stress (strain) displayed 
in Fig. [UJb). In the calculation shown here and after 
we take the pinning potential Ui to be set by the Debye 
temperature &d — 28.8 K [13, H|, since we postulate 
that pinning originates from dislocation crossing or local 
melting. Most materials parameters have been reported 
in the literature: b = 0.364 nm, ^ = 13.7 MPa [H], 
Tpeieris = 10~^ /i [20| . The fundamental results in this 
study are not sensitive to the assumed model parameters, 
such as a = 0.1, /? = 0.9 and c = 0.5. 

The slope of the stress-strain curve defines the work 
hardening rate (WHR= dTapp/dcapp). The WHR and 
dissipation in Fig. [3] show a decrease in the WHR above 
~ 100 mK accompanied by a dissipation peak. This re- 
sult is strikingly similar to measurements of the dynamic 
shear modulus. At low temperature the WHR is the 
value of an ideal crystal with no dissipation, since dislo- 
cation glide is negligible. At the crossover temperature, 
dislocations start to overcome barriers and glide. Dislo- 
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FIG. 4: (Color online). Work hardening rate versus temper- 
ature for different pinning point ('^He) concentrations. The 
curves for high, medium, low ^He concentrations correspond 
to initial mobile dislocation densities of 1 x 10*, 5 x 10*, 9 x 10* 
m~^, respectively. All calculations have the same initial to- 
tal dislocation density of 10® m~^ and total strain amplitude 
e = 1.7810 x lO"'^ 111. 



cation glide contributes additional strain, which reduces 
the WHR while simultaneously dissipating energy giving 
rise to a dissipation peak. At even higher temperatures, 
dislocation lines are essentially free to glide viscously; the 
system exhibits plastic strain but dissipates only mini- 
mal energy. We also show the WHR and dissipation for 
different applied strains in Fig. [3l The crossover temper- 
ature for a larger applied stress (strain) is lower since the 
"effective barrier" is lower in that case. 

Our dislocation glide model can capture the effect of 
■^He impurities as well. We propose that '^He atoms dec- 
orating dislocation lines exert a drag on moving disloca- 
tions, instead of pinning them since dissolved '^He atoms 
exhibit a very high mobility in solid "^He compared to 
the speed of dislocations. We model a lower ■^He concen- 
tration by a higher initial mobile density Pmob and vice 
versa for the same initial total dislocation density. Fig. |4] 
shows the WHR versus temperature. The curves for high, 
medium, low '^He concentrations correspond to initial 
mobile dislocation densities of 1 x 10^,5 x 10^,9 x 10^ 
m^^, respectively. When the '^He concentration is low, 
the main obstacle for dislocation gliding is still the Peierls 
potential and the background interaction. Hence the re- 
sult does not depend much on the "^He concentration. 
When the ■^He concentration is high enough, the viscous 
drag starts to dominate over other effects and causes a 
broader crossover at higher temperatures. This is in ac- 
cordance with the observations made in various experi- 
ments and shown for low, medium, and high '^He concen- 
trations in Fig. m 

We next discuss the similarities and differences be- 
tween our calculations and the previous AC strain mea- 
surements of the dynamic shear modulus [H-Q ■ Although 
we model an applied DC strain rate measurement, the 
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results share the same quahtative features with the AC 
experiment: both show decreases in the WHR above 
a crossover temperature accompanied by a dissipation 
peak. However, the magnitudes of these features be- 
tween our theory and measurements are quite different. 
Large differences in stress-strain relations have been re- 
ported before between free-standing and confined crystals 
of solid helium and may be a contributing factor to this 
discrepancy [2l|. The critical strain in the AC experi- 
ment is about 10~^, while in our model it is roughly 10~^. 
The ratio of change in the WHR is also much smaller in 
the AC measurement 15%) than in this calculation 
(~ 60%). This suggests that the AC experiment may 
measure a type of dislocation motion that dissipates less 
energy than gliding. In particular this motion should not 
involve moving the entire dislocation line from one Peierls 
valley to another, since the critical stress in the AC mea- 
surement (tc ~ 10~^/i) is at least two orders of magni- 
tude smaller than the Peierls stress {rpeieris > 10^^/i). 
A potential explanation for such motion includes dislo- 
cation vibration [l^ and kink motion [g^, H^- If this 
mechanism can be proven dominant, then the disloca- 
tion motion is different in the AC and DC strain rate 
experiments and may resolve the inconsistencies between 
the torsional oscillator (TO) results with and without a 
rotating DC component [13]. In this case, our dislocation 
glide model may be generalized to include dislocation vi- 
bration at low strains to explain the AC shear modulus 
and TO experiments with a superimposed DC rotation. 

Finally we note that grain boundaries and the dis- 
tribution of multiple grain orientations can change the 
strength of the effect. In a polycrystalline system, the 
strain from different grains can act along different di- 
rections and partially cancel each other out. This would 
diminish the observed effect. Indeed ultrasound measure- 
ments observed a smaller change in the shear modulus of 
polycrystalline samples compared to single crystals. For 
high quality single crystals a change of 86% was reported 



fS]. Our model does not include such cancellation ef- 
fects. This may partly explain why the AC strain rate 
measurement in polycrystals observes a smaller effect. 

In summary, we proposed a simplified isotropic 
density-based dislocation model to describe qualita- 
tively many essential features of the anomalous low- 
temperature shear modulus of solid "^He measured with 
an AC strain rate. More importantly, we predicted many 
new stress-strain curves, work hardening rates, and dissi- 
pation properties for a DC strain rate shear experiment. 
In our model, the classical glide of dislocations in slip 
planes under load has no relationship with supersolidity. 

One may speculate that the existence and dominance 
of dislocation induced plastic deformation facilitates a 
nonuniform supersolid state, because a supersolid state 
is possible in ''He along the core of static dislocations 
[25I [26| . An interesting aspect of a highly tangled dis- 
location network may be its local loss of crystalline or- 
der, due to strong lattice distortions by closely arranged 
dislocation cores [27|. Those distorted regions are more 
amenable for superfluidity and are a potential candidate 
for a nonuniform supersolid phase. The dynamic effects 
of plastic deformations and dislocation glide discussed 
here will play a role in the dynamic measurements of 
solid ''He regardless of the presence or absence of super- 
fluidity. In the absence of more detailed microstructural 
information the model proposed here offers an effective 
method to check whether the dislocation motion is the 
fundamental mechanism for plasticity in solid helium. 
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